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ABSTRACT 

A graphical representation of angular momentum is used to evaluate 
relativistic matrix elements between antisyrranetrized states of many- 
particle configurations having any number of open shells. The antisym- 
metrized matrix element is expanded as a .sum of semi symmetrized matrix 
elements, which can be evaluated expediently in terms of radial integrals 
from easily constructed diagrams. The diagram representing a semi sym- 
metrized matrix element is composed of four diagram blocks, namely, the 
bra block, the ket block, the spectator block, and the Interaction block. 
The first three blocks indicate the couplings of the two interacting con- 
figurations while the last depends on the interaction and is the re- 
placeable component. Interaction blocks for relativistic operators and 
commonly used .potentials are summarized in ready- to-use forms. A simple 
step-by-step procedure is prescribed generally for calculating anti sym- 
metrized matrix elements of one- and two-particle operators. A modified 
covariant 3-jm coefficient is also introduced along with certain new 
graphical notations. Although we focus on jj-coupled states, which comes 
naturally in a relativistic formulation, the general procedure holds in 


any coupling scheme. 
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I. INTRODUCTION 

A relativistic treatment Including quantum-electrodynamic corrections 
has become essential in many fields of atomic physics; this requires the 
evaluation of matrix elements with considerable complexity. In particular, 
the possibility of forming highly ionized heavy atoms and the study of 
their structure and of their collision and subsequent cascade processes 
necessitate a complete treatment of couplings between several open shells. 

On the other hand, the application of Wigner-Racah‘s idea in these problems 
•with standard techniques of angular-momentum coupling becomes a tedious and 
more often arduous task. The purpose of this work is to provide non- 
specialists with a simple and powerful tool to express complicated- matrix 
elements in terms of radial integrals, suitable for numerical computation. 

One of’ the essential approximations in the quantum-mechanical descrip- 
tion of a many-particle system is the central field approximation. Orbitals 
of the particles can thus be represented by angular-momentum eigenstates. 

The coupling of angular-momentum eigenstates with irreducible tensor 
operators depend only on the rotational -properties of the states and oper- 
ators involved. This fact leads naturally to the division of the calcula- 
tion of a physical quantity into, two parts: One consists of dynamical 

variables invariant under rotations, and the other is a geometrical factor 
depending on the rotational properties of the physical quantity. It is 
the Wigner-Eckart theorem (Wigner, 1927; Eckart, 1930) which embodies this 
notion. The geometrical factor is given by the Clebsch-Gordan coefficient 
(also called the Wigner or vector-coupl i ng coefficient ) or by a more sym- 
metric quantity, the Wigner 3-j coefficient (Wigner, 1951; Edmonds, 1957). 
Techniques for the solution of related algebraic problems for many-rparticle 
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systems have been developed by Racah {1942a, 1942b, 1943)’. [Two other 
equally applicable approaches have been described in the works of Condon 
and Shortley (1935) and Harter .and Patterson (1976), respectively, but 
will not be discussed here.] However, the complicated algebraic manipula- 
tions have prevented' nonspecialists from carrying out specific calcula- 
tions. 

Attempts to solve the complexity of angular-momentum couplings result ■ 
in graphical methods- In the early stage, angular-momentum d.iagrams were 
used only in a qualitative and descriptive way (see, e.g., Edmonds, 1957 
and Judd, 1963). It was the achievement of Jucys, Levinson, and Vanagas 
(1962) that a graphical method was put on a quantitative basis so that 
angular-momentum couplings can be solved in an expedient and elegant manner 
solely in terms of diagrams. The graphical method of Jucys et al . starts 
by assigning a graphical symbol to the Wigner 3-j coefficient and compounds 
the angular-momentum coupling to 3n-j symbols . An alternative method given 
by Danos (1971) however focuses on the recoupling aspect and uses the 9-j- 
recoupling as the central graphical element. The former has been particularly 
effective in extending the range of application of Racah 's techniques. The 
graphical approach permits transformations on diagrammatic expressions and 
leads on to analytical results in a clear and simple manner. Besides its 
utility as a calculational tool, the graphical method has the appealing 
feature of revealing, at a glance, the structure of very complicated 
couplings of angular momenta. In addition, we do gain in the graphical 
form some physical insight, similar to the visual understanding of physical 
processes provided by Feynman diagrams. In fact, because of the role which 
the graphical method plays in extracting the geometrical part of a Feynman 
diagram (Bolotin, Levinson, and Tolmacher, 1964; Judd, 1967; El-Baz and 
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■Nahabetlan, 1969), it becomes an indispensable supplement to the Feynman 
diagrams, where a perturbed quantum mechanical system is studied graphically. 

The graphical method has been developed subsequently by Jucys and 
Bandzaitis (1967), Hassot, El~Baz, and Lafoucriere (1967), Brink and 
Satchler (1968)., El-Baz (1969), Bordarier (1970), Briggs 0971), El-Baz 
and Castel (1971, 1972). El-Baz and co-workers have in particular extended 
the graphical method to treat spherical harmonics, irreducible tensor . 
operators, and rotation matrices. The basic idea of El-Baz and Castel (1972) 
consists in representing graphically the bra (covariant) and ket (contra- 
variant) vectors familiar in the Dirac notations C0i^‘ac, 1930). The 
graphical representation of the Clebsh-Gordan coefficient thus becomes a 
straightforward extension of the bra and ket diagrams. The graphical repre- 
sentation of the Wigner 3-j coefficient, first introduced fay Jucys et al . 
(1962), has been modified by El-Baz (1969) to better represent its co- 
'Variant property, using Wigner’s covariant notation (Wigner, 1959). We 
will adopt this idea of El-Baz with a modified phase f^tor to have a more 
coherent correspondence between the Wigner 3-j coefficient and the Clefasch- 
Gordan coefficient. 

A graphical treatment of antisynmetrization for the evaluation of anti- 
symmetrized matrix elements has been given fay Bordarier (19701 and Briggs 
(1971). The Bordarier’s treatment is general and encompasses many different 
types of matrix elements, while Briggs gives a step-by-step procedure for 
the evaluation of matrix elements between antisyimnetrized states with LS 
coupling'. An alternative treatment has been given by Huang and Starace 
(1978) for a particular case. The procedure of Bordarier and Briggs is 
however a little intricate between purely graphical steps and manipulations 
which are. better to be performed analytically. For exanple, the antisym- 

■r. — '^iffprent subshells may be carried out 
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analytically with ease without resorting to a graphical phase rule. Fur- 
thermore, the interaction diagram is obtained in the Briggs' prescription 
by expanding the interaction operator in a complete set of particle or- 
bitals, whereas a similar interaction diagram may be obtained by con- 
sidering directly the m-scheme matrix element. 

In this work, we will prescribe a simple step-by-step procedure for 
evaluating anti symmetrized matrix elements for one- and two-particle oper- 
ators, Although we focus on jj-coupled states, which comes naturally in 
a relativistic formulation, the general procedure holds in any coupling 
scheme. The underlying idea of the present approach is to express analyti- 
cally the matrix element betweisn antisymmetrized many-particle states in 
terms of matrix elements between semi symmetrized many-particle states . 

The semi symmetrized many-particle state Is defined as the many-particle 
state which is antisymmetric within each subshell but is not antisymmetric 
with respect to exchange of two particles from different sufashells. Hence- 
forth, the matrix elements between semi symmetrized many-particle states 
are evaluated by a graphical procedure. The graphical procedure consists 
of three major steps: Firstly, we construct diagrams for the two in- 
teracting semi symmetrized many-particle states and decouple active 
particles from the other particles to be referred to as the spectator 
particles . Here the active particles represent the typical particles 
which actually participate in the interaction in a semi symmetrized matrix 
element. Secondly, we insert the interaction block corresponding to the 
interaction between active particles. The interaction blocks for conmonly 
used operators and potentials are summarized in ready- to-use forms in 
Appendix B. Lastly, we evaluate the resultant diagram by transforming it 
into basic diagrams representing 3n-j symbols, whose analytical values have 
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been tabulated extensively {Rotenberg et al., 1959; Jucys et al., 1962). 
Thus the anti symmetrized matrix element is expressed as a sum of weighted 
radial integrals. 

The graphical notation and transformation rules' used in this worR are 
given in Sec. II. The covaria'nt 3-jm coefficient is defined there. In 
Sec. Ill an analytical procedure is outlined for evaluating the antisym- 
metrized matrix element in terms of radial integrals. In Sec. IV we 
describe in detail how the diagram representing the semisymmetrized matrix' 
element can be constructed and, eval uated. In Sec. V we summarize the pro- 
cedure of evaluating antisymmetrized matrix elements. A worked example ts 
given in Sec. VI. Appendix. A is a glossary of definitions of terms used 
in this' work; in most cases they are identified in the text by italic 
letters when they first appear. Graphical forms of commonly used operators 
and potentials are surmiarized in Appendix B. Appendix C gives the deriva- 
tion of one of the expansion formulas used in Appendix B. 
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II. GRAPHICAL NOTATION AND TRANSFORMATION RULES 

We will first define the basic coupling coefficient, i.e,, the 
covariant 3-jm coefficient, and its graphical representation in Sec. 11. A. 
Rules to combine 3-jin diagrams, which amount to sunwations over magnetic 
quantum numbers m, and rules to transform the combined diagram will be 
given in Sec. II. B. Section II. C suirmarizes analytical values of some 
basic diagrams. The graphical representation of the Cl ebsch-Gordan co- 
efficient and its transformation to the 3-jm diagram is described in ■ 

Sec. II. D. These find uses In constructing configuration diagrams. 

A. The Covariant 3-jm Coefficient 


Wigner's 3-j symbol is defined by 



X 


- 




= C-) 




1.2.n 


where <jimi j 2 ni 2 l is the Clebsch-Gordan coefficient. We employ, however, 
the covariant .notation (Wigner, 1959) such that the covartant 3-jm coefficient 
(or, simply, the 3-jm coefficient) which is covariant in the first two in- 
dices and contravariant in the last index is defined- by 


i K T » 

Wi, 3^ J 












( 2 . 2 ). 



10 


Note that this definition is different from Wigner‘s definition CWigner, 

j 

1959) in that he used the phase factor C“) • This is also different 

from the definition of El'-Baz and Castel (1972) in that our contravariant 
component is defined with the same phase as their covariant component. 
This modification, however, does not change the graphical rules presented 
in their work. We also emphasize that th.is coefficient is called "the 
covariant 3-jm coefficient" because of their covariance properties and 
m-dependence. Accordingly the name "3-j symbol" or "3-j coefficient" 
will be reserved for the coefficient occurring in the hierarchy of 3n-j 
coefficients, which have no m-dependence. 

By our definition, the ordinary 3-j coefficient of Wigner is a fully 
covariant 3-jm coefficient and is equivalent to the fully contravariant 



(2.31 


C2.4) 


The notational rules are as follows: 



(i) Each vertex, indicated by a node, represents a 3-jm coefficient. 
Each covariant component is denoted fay a line with an ingoing double-arrow, 
and each contravariant component by a line with an outgoing single-arrow. 
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(ii) The plus (minus) sign at the node means that the angular momenta 
are to be read counterclockwise (clockwise). The change of the sign at the 
node (3102^3) introduces a phase factor (-) ; therefore, the sign at 
the node may be suppressed if number. 

The magnetic quantum number of each angular momentum line is usually 


suppressed where no confusion may occur. Also, we note that the "current" 
of magnetic quantum numbers is conserved at each node due to the selection 
rule for magnetic quantum numbers of the covariant 3 -jm coefficient; for 



where we have the sum of the ingoing currents "m-j" + "mg" = the outgoing 
current "m-| + mg." 


The notational rules used to combine 3 -jm diagrams are as follows: 

(i) The summation, or in the tensorial terms "contraction," over a 
pair of magnetic quantum numbers (one of which is always contravariant, and 
the other covariant) is performed by joining the corresponding angular- 
momentum lines to form a linked single-arrowed line. ' 

(ii) The change in direction of a linked angular-momentum line j 

21 

introduces a phase factor (-) As a result, we may suppress the arrow 
of a linked angular-momentum line j whenever. j is an integer. 

Although we do not write the magnetic .quantum number for a Tinked 
angular-momentum line, the summation over the magnetic quantum number is 
always implied. Nevertheless, in many cases because of the conservation 
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of the "magnetic current," the summation implied by a linked angular- 
momentum line exists only formally. - 

B. Transformation Rules for 3-Jm Diagrams 

There are only two fundamental transformation rules: 


Rule I: 


where a bar on the angular momentum k stands for a multipTi cation factor 
(2k+l)^, and multiple bars for multiple .factors. For example, R indicates 
a factor (2k+l). Rule I follows from the graphical relation for 3-jm 
coefficients. 





which represents the orthogonality relation 



In C2.7) we use the blocks a and 6 to represent arbitrary diagrams either 
open or closed. By open or closed , we mean that the diagram either has or 
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hasn't any free nonzero angular-momentum lines. For example, In the fol- 
lowing diagram 



the diagram block on the left Is closed while the one on the right open. 
We will use an encircled greek letter © to indicate specifically a 
closed diagram block. 

Rule II: 




C2.10) 


where the bar under the .angular momentum j, denotes a multiplication factor 
(20i+l) . Also we will use multiple bars to' indicate multiple factors. 

Since derivation of this transformation rule is more Involved, we refer the 
reader to the works of Jucys et al. (1962) and El-Baz 'and -Castel (1972) for 
its proof. Note that for a null block e. Rule II becomes 


© 

■ — ^ — 




0 






■ 






( 2 . 11 ) 
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This transformation rule results from the rotational invariance, of the 

¥■ 

diagram. 

From these two fundamental rules, we can easily derive besides others 
the. following additional useful rules: 
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The transformation rules I and li and the rules derived from them 
allow one to rearrange diagrams or to factor out basic diagrams whose 
analytical values have been tabulated. Some of those basic diagrams are 
presented in the next subsection. 


C. Analytical Values of Some Basic Diagrams 




C2.161 


C2.17) 



( 2 . 20 ) 
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+ 





In general, we have 



which is defined as the Sn-^j symbol of the first kind . For-n = 1, 2, 3, it . 
equals the ordinary 3n-j symbol within a phase factor. We also have the 
3n-j symbol of the second kind 



(2'.23) 
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Further discussion of these 3n-j symbols may be found in the works by Jucys 
et ale (1962) and by El-Baz and Caste! (1972). We emphasize that it is the 
topology of an angular-momentum diagram how the various lines in it 

are connected to each other) .which determines its analytical value. There- 
fore by keeping the topology, we can deform an angular-momentum diagram in 
any way without changing its analytical value. 

A very complete tabulation of Clebsch-Gordan coefficients can be found • 
in Tables of the Clebsch-Gordan Coefficients , compiled, by the Institute of 
Atomic Energy, Academica Sihica (1965). Extensive tabulation of the 3-jm 
symbol and the 6-j symbol may be found in the work by Rotenberg et al. 

(1959). This reference also contains extensive references in the literature 
on 3n-j symbols. 

D. The Clebsch-Gordan Coefficient 

In constructing coupled angular-momentum states, where the Clebsch- 
Gordan coefficient (to be referred to as the C-G coefficient) occurs" 
naturally, one finds that the graphical representation of the C-G coeffi- 
cient is more expedient to- use. A graphical representation (El-Baz, 1969) 
of the C-G coefficient is given by 

The notational rules to represent and to combine C-G coefficients are the 
same as those for 3-jm coefficients. The only difference is that we use 
an open circle (instead of a solid- node as in the case of a 3-jm coefficient) 





12.24) 
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at the vertex to represent a C-G coefficient. Although transformation 
rules for C-G diagrams have been given (El Baz, 1969; El-Baz and Castel, 
1971* 1972)s we will not present them here because we only use the graphical 
representation of the C-G coefficient in constructing many-particle con- 
figurations. The resulting C-G diagram will then be transformed into a 
3-jm diagram by a. simple procedure. 

To transform a C-6 diagram into a 3-jm diagram, we do the following: 

(i) Add a bar, J, on the angular momentum j which is unique in its 
variance character. The bar represents a multiplication factor (2j+l)'^. 

(ii) The second (third) angular momentum k, counting from J in the 

direction indicated by the sign of the open circle, introduces a phase 

factor (-) if j is contravariant (covariant). We write a circle around 

the arrow of the angular-momentum line k to indicate this phase factor. 

When the angular-momentum line k is joined to another angular-momentum line- 

2k 

to' form a closed line, the circle representing the phase factor C-) . 
simply changes the direction of the closed angular-momentum line k. 

(iii) Fill the open circle at the vertex into a solid node and change 
the sign at the node. 

Note the special case 




with the sign at the node unchanged 


Similarly we have the equivalence 
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• We will find frequent uses of these simple relations, in constructing many- 
particle configurations. 
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III. ANALYTICAL EXPANSION OF ANTISYMMETRIZED MATRIX ELEMENTS 


The method of evaluating matrix elements between anti syrmetri zed many- 
particle states is outlined analytically in this section. The part of the 
manipulation which can be performed more expediently by a graphical proce- 
dure is presented again in detail with graphical representation in the next 
section. 

In the Dirac-Fock or Dirac-Fock-Slater description of a many-particle 
system, a configuration is specified by the number of equivalent particles 
occupying each subshell. In general, there are particles in the suBshell 

a 

jg, N|j particles in the sufashell . . . , particles in the subshelT 

j , etc. Accordingly this configuration is denoted by 
A 



For reference purpose, we may order the subshells in a certain sequence 
(e.g,, Is^, 2si^, 2p^, etc. for atomic subshells) and refer then by suc- 
cessive values of the index X. 

We consider states of the many-particle system in which each subshell 
of equivalent particles is in a definite state with total subshell angular 
momenta J,, . . . J, , etc. Within a subshell j. , if ‘ there are several 

QUA A • 

states with the same total subshell angular momentum there will be an 
additional quantum number which is required to specify a subshell state 
uniquely. This coupling scheme is indicated fay 






r 



All the total subshell angular momenta (J^ . . . ) are then 

coupled successively to form a grand total angular momentum J of the whole 
many-particle system; After anti symmetrizat ion, such a state is denoted by 

(3.1) 

where a stands symbolically for the coupling scheme of total subshell an- 
gular momenta. 

We will first express the matrix element. between anti symmetrized many- 
particle states (3.1) as a sum of matrix elements between semi synmetri zed 
many-particle states . As mentioned in Sec.. I, the semi symmetrized many- 
particle state is defined as the many-particle state which is anti synmetri c 
within each subshell but is not antisymmetric with respect to exchange of 
two particles from different subshells. Each semi synmetri zed matrix element 
is then expanded in terms of jm-scheme matrix elements . After the evalua- 
tion of jm-scheme matrix elements, we thus obtain the expansion of the 
anti symmetrized matrix element in terms of reduced matrix elements,, i.e., 
weighted radial integrals. 

A. Expansion of Anti symmetrized' Matrix Elements in Terms of Semi syranetrt zed 
Matrix Elements 

We consider matrix elements of one-particle operators and two-particle 
operators in turn. 

One-particle Operator 
We define the- one-particle operator by 




22 


From the general consideration of de-Shalit and Talmi 0963) and of Fano 
(1965), we can easily deduce the result 

|[( i_‘^) ^ 

= L <V"’> . . ■ 


i = i 


ft-b 


(3.3) 


where the summation is over all nonvanishing subshell pairs (a» h) with 
each pair counted once, and is defined as 


(3.4) 


with 


T -L N 

A-a.+ i ^ 


C3.51 


Here <q(a)aOH[ and jq‘ (h)a'J'H'> are semisymmetrized many-particle states 
with definite particle distributions specified by q(a) and q‘(b), respec- 
tively. Explicitly q(a) and q‘(h) denote the particle distributions in 
which the Nth particle (i.e., the active particle ) is in subshells a and 
b, respectively, while all the other particles (i.e., the spectator 
particles ) assume the same distribution in both states. We emphasize 
that how the spectator particles are distributed among subshells is 
immaterial so Tong as they keep the same distribution in both states. 
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Two-particle Operator 


We define the two-particle operator by 

I T/ . 

L<i 


C3.6) 


Also from the general consideration of de-Shalit and Talmi (1963.) and of 
Fano (1965), we can deduce the result 

^ j , » / ^ O if 

4- • •] 


E <v'‘> . , 

Ab, C ci * <Xb , t d. 


(3.7) 


where the summation is over all distinct nonvanishing pairs with a <_ b and 
( 2 ) 

c d, and <V' '>^ 1 ^ is defined as 


<V'"> 


= (-) 


fll' 


Abc4 


with 




C3.8) 


Lc.. = i 

c^kcd. >.= A+i ^ Ok=C+l 




C3.aj 
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Here <qUb)aJMr denotes a semi synsnetri zed many-particle state. with the 
(N-l)th and Nth particles (the active particles in subshells a and 
b, respectively. The distribution q'(cd) and q'(dc) are defined similarly 
to q(ab). Again we emphasize that all of them have the same spectator- 
particle distribution. 

B, Expansion of Semi symmetrized Matrix Elements in Terms of jm-Scheme 

Matrix El ements 

To evaluate matrix elements between semi symmetrized many-particle 
states we need to single out those particles which, actually participate 
in the interaction, i.e., the Nth particle and the CN-l)th and Nth par- 
ticles in (2.4) and (2.8), respectively. This may be accomplished 5y 
fractional parentage expansions of the subshell states involving active 
particles. With coefficients of fractional parentage (to be referred to 
as c.f.p.) as expansion coefficients, the semi symmetrized many-particle 
state can thereby be expressed by a linear combination of parent states. 
Each of these parent states can then be decoupled into a product of two 
parts: One contains active particles, and the other contains spectator 
particles. These expansions and decouplings enable us to express a 
semi symmetrized matrix element in terms of jm-scheme matrix elements. 

We consider the cases for one-particle operators and for two-particle 
operators in turn: 

One-particle operator 

To evaluate the semi symmetrized matrix element in (3.4) we first 
decouple the semi symmetrized many-particle states as 

r 


(3,10) 
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and 


|?(b)=(j'M')“ L C|,(p'; i'ib)} I t > . . 

f' 


C3.n.) 


Here <p;q(a)[ and |p';q'(b)> denote symbolically ”tne uncoupled subshell 
states of the spectator particles, i.e.,^the first (N-1 ) particles; <aj 
and |b>, or explicitly and are the orbitals of the active 

particle, i.e., the Nth particle. The expansion coefficient C,(p;aJM) or 

a 

C|jCp* ; a' J'M*') stands symbolically for the product of a c.f.p. and all the. 

3- jm coefficients needed in the uncoupling, and the summation index p or 
p‘ for- the summation over the c.f.p. and magnetic quantum numbers. Explicit 
examples will be given in Sec. IV. B when we consider the graphical proce- 
dure. By using {3.10)- and (3.11), we can write the semi synmetri zed matrix 
element in (3.4) as 

^=< T M 1 1/ I ^ V b ) J M ' > 




(3.12) 


Here the matrix element <p;q(a)| p' ;q‘ (b)> represents a product of overlap 
integrals and is independent of the interaction. The matrix element, 
<a|vj^|b>, called the jm-scheme matrix element , depends on the interaction 
and will be evaluated in the next subsection. Interested readers may refer 
to (4.15) for the graphical representation of (3.12). 
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Two-particle operator 


We decouple the semisynmetrized many particle states in (3.8) and 


P 


T' C3.14) 


Here <p;q{ab)[ and <p';q‘(cd)| stand symbolically for. the uncoupled subshell 
states of the spectator particles, i.e., the first (N-2) particles; 

-<ab| and [cd>, or explicitly <j^m^ and for the active 

particles, i.e., the (N-l)th and Nth particles; Cgjj(p;«JM) and C^^(p' ;a' J 'M' ) 
for the expansion coefficients which are products of c.f.p. and 3-jm coeffi- 
cients; p and p* for all the summation indices involved. Explicit examples, 
will be given in Sec. IV. B. Note that here a and b (also c and d) may 
represent either equivalent or nonequivalent orbitals. • 

By means of 'the expansions (3.13) and (3.14) we- obtain the first semi- 
symmetrized matrix element. (the direct term) in (3.8) as 

» (erf)) 

r1“ " / 

|cc(\ (3.15) 
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Here the matrix element <p;q(ab) j p* ;q'Ccd}> represents a product of overlap 
integrals, and <abj f|icd> is a jm-scheme matrix element. The second 
semisyrmetri zed matrix element (the exchange term) in (3.8) can be expanded 
similarly as (3.15) in terms of jm-scheme matrix elements. Interested 
readers may refer- to (4.16) for -the graphical representation of (3,15). 

C. jm-Scheme Matrix Elements 

In the last subsection, we have shown how a semi symmetrized matrix 
element can be expended as a sum of products of two parts: One is the 

interaction-independent part involving the coupling coefficients, and the 
other is the interaction-dependent part represented by a jm-scheme 
matrix element. The jm-scheme matrix elements for the cases of one- 
particle operators and two-particle operators are g,iven by <a|Vj^lb> 
and <ab|Vj^_.^ f^!cd>, respectively. For specific operators, these matrix 
elements can be evaluated analytically in terms of radial integrals. The 
results for commonly used operators .and potentials are presented in 
Appendix B along with their graphical forms. 

In general, any operator can be written as a sum of products of 
irreducible tensor operators. Hence in this subsection for. a general 
purpose, we consider Vjq and to be irreducible tensor operators. 

The results are given in terms of reduced matrix elements as follows. 

One-particle operator 

Assume the one-particle operator to be an irreducible tensor 
operator of degree j, 

'n = TjnW) . (3.16) 



By applying the Wigner-Eckart theorem (Wigner, 1927; Eckart, 1930) we 
obtain the jm-scheme matrix element in (3.12) as 


a- ,^‘b 

Here we denote the angular momentum coupling by a covariant 3-jm coefficient 
‘'jail 1 1 reduced matrix element which is usually expressibl 
as a sum of weighted radial integrals for a. specific case. It is of 
interest to note that the bra Ccovariant) state <j' m,i corresponds to the ■ 

a a ' 

covariant component in the 3-jm coefficient, and the contravariant operator 
^jm (contravariant) state ijfjm|^> correspond to the two contra- 

variant components. No extra phase or weight factor, besides the reduced 
matrix element, is carried by (3J7),and the rotational properties of the 
matrix element is clearly indicated by the 3-jm coefficient. 

To emphasize the fact that the jm-scheme matrix element is separated 
into a geometric part and a dynamical part, we rewrite C3.17) as 


>. 


(3.17) 



The geometric factor G.(a;b), which corresponds to a coupling diagram in . 

J 

the graphical representation and will be called the interaction graph , is 
defined in this simple case as 




The interaction strength X .(a;b) is defined accordingly as 

J 


C3.19) 
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Two-particle operator 

For the case of a two-particle operator, we. consider an irreducible 
tensor of degree j which is the tensorial product of two irreducible tensor 
operators acting on two different particles, i.e.. 


-If 

W-l, A/ 


(3.21)' 


As before, by applying the Wigner-Eckart theorem we can obtain the jm-scheme 
matrix element in (3.15) as 




= &• c.^) X- • 

tf* rf/ ^ 2. 


(3.22) 


Here the geometric factor G.(ab;cd) is given by 

vi 



iz y 


Mot 

^a. 


> 


(3.23) 


and the interaction strength by 



(3.24) 


D. Expansion of Anti symmetrized Matrix Elements in Terms of Reduced 
Matrix Elements 

- ' -In Sec, 1 1 1. A we expanded the matrix element between anti symmetrized 
“many-particle states as a sum of semi syiraietri zed matrix elements. In Sec. 
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III.B we obtained the semi symmetrized matrix element in terms of jm-scheme 
matrix elements, which were later evaluated’ in Sec. III.C. Here we summarize 
the results. 


One-particle operator 


For V' ^ -jai anti symmetrized matrix element (3.3) has 


been evaluated as 


where the coupling coefficient Dj„,Ca;b) is given by 


? 




(3.25) 


ft' 


X C‘<- ; l> J . 


(3.26) 


Two-particle operator 


,( 2 ) . " 


For = 2 T. (ij) with T. (ij) defined in (3.21),- we can surraiarize 

*j ^ j jin jrn 

the results for the anti symmetrized matrix element (3.7) as 

aS,c4. ^ 


• C3.27). 
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Here the coupling coefficient D. (ab;cd) is defined as 


L C^Cf J £Ktm)C . Cf f '(ccO 
tf' 

X Gf • ^ ^ 


with D. (abide) similarly defined.. The- interaction strength X. . (abicd) 

jHl '*1^2 

has been given in (3.24), and X. . (abide) is given by a similar expression 

^ 1^2 

with c and d interchanged. 
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IV. GRAPHICAL EVALUATION OF SEMISYMHETRIZED MATRIX ELEMENTS 

Evaluation of the coupling coefficients and Dj.j^Cabjcd) , given 

in (3.26) and {3.28), respectively, is a formidable task. Although for 
many cases they may be evaluated numerically CG>^ant, 1973; 1976)' by using 
a digital computer, these coupling coefficients in general have to be ob- 
tained analytically for each particular case, especially for the analytical 
Study of a matrix element. In this section we will show how to obtain an 
analytical expression of an anti symmetrized many-particle matrix element 
from easily constructed diagrams. 

In Sec. IV. A the graphical procedure of constructing semi symnetri zed 
many-particle states is illustrated by an example. Section IV. B describes 
how to decouple graphically particles from an anttsynmetrized subshell 
state. The bra and ket diagram blocks are also defined there. In Sec. IV. C 
the jm-scheme matrix element is considered. Its graphical representation is 
defined as the interaction block . Specific diagrams are given for the cases 
considered analytically in Sec. TII.C. Other interaction diagrams for 
commonly used operators and potentials are presented in Appendix B. In 
Sec. IV, D the spectator block is defined. Evaluation of the joined diagram, 
called the recoupling diagram , is described in Sec. IV.E.‘ 

A. Construction of Semi symmetri zed Many-Particle States 

The construction of diagrams for semi symmetri zed many-particle states cai 
easily be carried out in the C-G representation. The procedure is best 
demonstrated by working with an example. 

Consider a configuration having open shells a, b, c, d, and closed 
shells A, etc. . .. A particular coupling scheme of the open shells is rep- 
resented by 
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|[<Vb)Jab'JcJd)JcdW’ 


where J^, J^, are the total angular momenta of respective subshells, 

« 

and the parentheses specify the sequence of the couplings. These couplings 
can be given analytically by 


IC<Vb»ab(Vd)'>cdW> 

ltVb)Jab«abHCVdWcd«cd>'JabVcd«cdl*^ 

”ab"cd 

*^'>a»a'>bNblJab»ab=-<''cVdMdNcd«cd><'^ab«abJcd«cdlJ»> • W-D 


which corresponds to the graphical representation 














4 


C4.2a) 
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or simply. 


C4.2bl 

All the closed shells are represented symbolically by 

T^ = 0- 

I 



Hence the semi synimetri zed many-particl e state has the ciraphical represen- 


tation 



C4.3I 
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where a denotes symbolically the coupling scheme of the open shells. 

From the above example we can generalize the following graphical pro- 
cedure for constructing the semi symmetrized bra and ket states: 

(i) Represent the grand total angular momentum by- a line. 

(ii) Draw two branches from the grand total angular momentum and 

a small circle at the vertex. The two angular-momentum lines added repre- 
sent the two angular-momenta which are coupled to form the grand total an- 
gular momentum. Write a sign, + or at the vertex to indicate the order 
of coupling. 

(iii) Repeat step Cii)» starting from the new angular-momentum lines, 
until all subshells containing active particles are decoupled. 

(iv) Add angular-momentum lines for closed shells. 

(v) Mark appropriate arrows on the angular-momentum lines. 

It is more convenient -to construct bra (covariant) states from left to 
right and ket (contravariant) states from right to left. The so-constructed 
configuration diagrams are to be transformed into diagrams in the 3-jm 
representation by using the procedure given in Sec. II. D. 

B, Decoupling of Active Particles; Bra Block and Ket Block 

As stated in Sec. III.B, we can single out active particles from a 
subshell by a fractional parentage expansion. We consider two cases in 
turn. 


One-particle coefficient of fractional parentage 

To single out one particle from a subshell we use a one-particle frac- 
tional parentage expansion, i.e.. 
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C‘^j"“->'X.T,3-('‘T) <T,M, ) 

«^,T, 

Here to indicate more clearly the coupling, we use the abbreviated 
notation 

for the c.f.p. defined by Racah 0943).. 

The last expression in (4.4) gives the explicit form “of '(3.11) in the 
particular case of one subshell. This decoupling can be represented fay • 
the diagram 



C4.6) 


^et '^ocA. 
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where the c.f.p. [{j'^"bct^J-j j |aJ] and the sumnation over it are Implied by 
the square [U at the vertex (J-| j J). The ket block is defined as the por- 
tion of the diagram indicated in (4.6); a bra block' is similarly defined in 
the case of a covariant state. 

Two -particle coefficient of fractional parentage 


To single out two particles from a subshell we use a two-particle 


fractional parentage expansion, i.e.. 


■>n w' 







C4.7). 


where we have used the abbreviated notation 

' (4.8) 



to represent the two-particle c.f.p. The last expression of (4.7) gives 
the explicit example for separating two equivalent electrons from a sub- 
:Shell, which was given s^bolically in (3.14). 

The graphical representation of (4,7) is given by 


I ^ 



(4r iXoJk. 


C4.9) 


where the two-rparticle c.f.p. [(j” )a2«32Cj summation over 

it are denoted by the double squares at the vertex (JgJ’J). The ket block 
is defined as indicated in (4.9); a bra block is defined in the case of 
a covariant state. 

Tables of c.f.p. may be found in the works of Edmonds and Flowers (1952) 
de Shalit and Talmi (1963), and Sivcev et al . (1974). . 

C. Interaction Block 

An interaction block refers to the diagram block representing <a|vj^l&> i 
the case of a one-particle operator and <ab|Vj^_^ {^lcd> in the case of a two- 
particle operator. These jm-scheme matrix elements have been given analyti- 
cally for general cases in Sec, III.C. Besides a dynamical multiplication 
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factor, the interaction strength , the jm-scheme matrix elements depend 
only on the rotational properties of the interaction and states involved. 
We present here the interaction blocks for the general operators worked 
out in Sec, III.C, Specific examples will be given in Appendix B. Again 
we consider two cases in turn: 


One-particle operator 


An elementary interaction is represented by an irreducible tensor 
operator of degree j. Its jm-scheme matrix element was given in (3.18) 
and has the graphical representation 


I > = 



(4.10) 


where the coupling coefficient G.(a;b.) is represented by the coupling diagran 

wl 

and the interaction strength X.(a;b) is denoted by the cross "X" at the verte 
(j’a 3 jjj). More complicated Interactions can be expressed as a linear com- 
bination of this elementary interaction. 

Two-particle operator 

We assume the elementary two-particle interaction to be the irreducible 
tensor operator N) defined in (3.21). Hence from (3.22) we obtain 

the graphical representation 
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C4.ll) 


Here the cross "X" at the vertex Cj-j j ^2^ denotes the interaction strength 
Xjij^(afa;cd); the bar on the angular momentum J represents the multiplication 

factor (2j+T)^. 


D. Spectator Block 

A spectator block refers to the diagram block representing the scalar 
product of subshell states or groups of subshell states which do, not partici- 
pate In the interaction considered in the semi synmetri zed matrix element. It 
examples were given as <p;q{a)lp' ;q‘ (b)> in C3.12) and <p;qCab) jp' ;q*Ccd)> 
in (3.15). As mentioned in Sec. III. A* the spectator particles, have the 
same distribution among subshells in both the bra and ket states. This fact 
Implies that the scalar product of the composite states can be written as a 
product of overlap integrals for all subshells. Graphically, we represent 
spectator blocks as follows. 



41 



(4.12) 


where and parent states of subshells ^ and ^ after the 

active particle being decoupled out;-a^J^ and denote symbolically all 
the other subshell states or groups of s'ubshell states. From (2.16) we know 
that (4.12) represent a product of Kronecker deltas for the states involved, 
provided we use the same orthonormal set of particle orbitals in both the - 
bra and ket states. - 


(11) <p;q(ab)|p';q'(cd)> 



( 4 . 13 ) 
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where the notations are defined similarly as in the case (1). Another 
example is given by 

f'; 



C4.14) 


where «j,2'^a2 state of subshell Rafter two active particles being 

decoupled out by a two-particle fractional parentage expansion. 

The graphical rule implied by (4,12), (4.13), and (4.14) is that we 
simply join together the corresponding (covariant-contravariant) angular- 


momentum lines in the bra and ket blocks. Note however that this -simple 
graphical rule does not apply when particle orbitals in the bra state are 
not orthonormal to those in the ket state. In such cases we need an extra • 


factor which is the product of all the overlap integrals of the spectator 
particles of the bra and ket states. 
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E. Evaluation of Recoupling Diagrams 

In Sec. IV.A we demonstrated how to construct diagrams for semisynmetrizi 
many-parti cl e states, and a simple graphical procedure was given. In Sec. 

IV. B we showed how to separate graphically active particles from a subshell. 
There, the part of the configuration diagram involving expansion coefficients 
and angular-momentum couplings constituted the bra block or the ket block in 
the case of a covariant state or contra variant state,' respectively. In Sec. 
IV. C and IV. D, we defined the interaction block and the spectator block; 
their typical diagrams were given. Here we summarize the results by the 
symbolic diagrams: ... 

(i) <q(a)aJM[v^jq'(fa)a'J‘M‘> 


^pccfbfbr' 





(. 4 . 15 ) 
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(i1) <q(ab)aJMiV[^|_-,^jjiq'(cd}a*J'M'> = 



C4J6). 
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After extracting the Interaction strength, c.f.p., and summations- over them i 
(4.15) or (4.16), we are in general left with a pure angular-momentum- 
coupling diagram, i.e.. 



. C4.17) 


Here the double-lines stand symbolically for all angular-momentum lines 
connecting two diagram blocks. By using the transformation rule (2.14), we 
obtain 



(4.18) 


The first factor in (4.18) represents a 3-jm coefficient, which is to be 
expected by applying directly the Wigner-Eckart theorem to the semi symmetri zee 
matrix element. The second factor in (4.18) is a recoupling diagram repre- 
senting analytically a recoupling coefficient. By using transformation 
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rules given In Sec. II. B, we can express the recoupling diagram in terms 
of products of 3n-j diagrams whose analytical values have Been tabulated 
extensively (Rotenberg et al., 1959; Jucys et al., 1962). 

The most expedient way to factor graphically a recoupling diagram 
into 3n-j diagrams depends, of course, on the particular diagram in 
question- A simple rule is to look first for diagram blocks separable on 
one angular-momentum line and then on two and three angular-momentum lines. 
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Ve SUMMARIZED PROCEDURE FOR EVALUATING ANTISYMMETRIZED MATRIX ELEMENTS 

In this section the results of Sec. Ill and IV will be suinmarized in 
the form of a step-by-step procedure by. which matrix elements of operators 
or potentials between anti symmetrized many-particle states can be evaluated. 
The prescription is given as follows. 

(i) Follow (3.3) and (3.4) for one-particle operators, or (3.7) and 
(3.8) for two-particle operators to express the anti symmetrized matrix 
element as a sum of semi symmetrized matrix elements. 

(ii) For each semi syiranetri zed matrix element, construct graphically 
semi symmetrized many-particle states; 

(a) Begin with the grand total angular momentum and the two 
angular momenta coupled to it. Here we use the bra state 
to illustrate the procedure: 



C5.1) 

(b) Repeat step Cal starting from and until all active 
suhshells are decoupled. Closed shells are then added 
separately. The resultant diagram is giVert schematically 
as 
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(.5.2) 

Note that in constructing semisymmetrlzed many-particle states those sub- 
shells which neither involve coupling with other subshells nor contain 
active particles may be ignored. 

The ket state is- constructed similarly. 

(iii) Decouple active particles from respective subshells by making 
use of the fractional parentage expansions (4.6) or (4.9), e.g.. 


(5.3) 

(iv) Repeat steps (ii) and (iii) for the ket state- 

(v) Insert the interaction block between lines representing active 
particles, and connect the corresponding spectator-particle lines in the bra 
and ket states. Here the interaction strength in the interaction block can 
be obtained analytically. A summary of the interaction strengths for various 
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conmonly used operators and potentials is given in Appendix B. 

(vi) Transform the C-6 representation of the bra and ket states into 
the 3-jm representation by using the transformation procedure presented in 
Sec. II. D. 


(vii) Transcribe analytical expressions from symbols representing the 
interaction strength, the c.f.p.'s and the summations association with them. 
The rest of the diagram represents a standard recoupling diagram which may 
be reduced to an analytical expression. The transformation rules to reduce 
an arbitrary recoupling diagram to tabulated 3n-j coefficients is given in 
Sec. II. B. 


(viii) Repeat steps Cii)“(vii) for other non-vanishing sufashell pairs, 
as indicated in step (i). The final expression of the anti symmetrized matrix 
element is thus obtained by summarizing all the analytical expressions re^ 
suited in step (vii). 
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VI. EXAMPLE 

We illustrate the graphical procedure by an example: 

Consider the C-AB;D Auger transition (Auger, 1925; Bambynek et al., 1972 
in a rare gas ion, where C denotes the initial vacancy in subshell C, and the 
final state is characterized by one vacancy in each of subshells A and B 
plus an outgoing electron D. This Is one of the de-excitation- processes of 
an atom with an inner-shell vacancy, in which the transition energy from 
filling the inner-shell vacancy by an outer-shell electron is carried off by 
the ejection of another outer-shell electron. The transition probability am- 
plitude is given in the Dirac-Fock formulation by (Huang, 1978a) 

- A6jP) ='([ <■ j-c ■=' I 

x-iv. 

: '-i 


C6-.ll 


’ ( 6.21 

where the bra state is the initial state and the ket state is the. final state 


with 


V.. 


/ - ^ e 

-7 -- «j) — 
• ' 'J 0 


t It? r. . 
•'J 




+ (.«; -V.K-X.-Vj) — 


’J-l 
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with the angular momentum of the continuum electron. Here all the other 
closed shells are denoted symbolically by <3^. To focus on the essential 
feature, we assume that the same orthonormal set of single-particle orbitals 
is used for both the final and initial states. By applying the step-by-step 
procedure in Sec. V, we evaluate the anti symmetrized matrix element (6.1) as 
follows. 


Step (i): There is only one nonvanishing subshell pair, i.e.. 


t(c^a0, .!>) = <£ \/^. > 


with 


l.<; 




C6.31 




Nb = 4- 1 






^c ’ 




N- = 2j^ + 1 


Nj = l 


(6.4) 


and 


? == 

'^A.bcoL 




I + ^ Gve»t 


)-a 
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Hence from (3.8) we can rewrite (6.3) as 

T^Cc-a$.5)= 

C6.6). 

The coupling schemes, indicated implicitly 5y a and a' for the initial 
and final states, respectively, have not been specified yet. The initial 
state is composed of closed shells with a vacancy in one of them and there- 
fore has no term structure. The final state consists of two vacancies and 
one continuum electron. We consider the coupling scheme such that the two 
almost filled subshells A and B are first coupled to form an ion core speci- 
fied by the total angular momentum The core state is then coupled to 
the continuum orbital D. This coupling scheme of the final state is given 
explicitly as 

Now we proceed to evaluate the semi symmetrized matrix elements in (6.6) with 
the coupling scheme (6.7) for the final state. 



53 


Direct matrix- element 

For the first term in the curly brackets in (6.6), we obtain the 
following. 


Step (ii)-(v): Because those subshells which neither involve 

coupling with other subshells nor contain active particles can be ignored, 
we obtain 



IrCjU j. atatX 


IrttiWActt SIX 
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Here the interaction strength, denoted by the cross "X" in the interaction 
block, of the potential (6.2) can easily be looked up in Appendix B. Note 
that all the c.f.p. in the diagram (6,8) are trivially unity although we 
retain the c.f.p. symbols to illustrate the decouplings. Also we have made 
use of the simple relations (2.25) and (2.26) in decoupling closed shells. 

Step (vi): The 3-jm representation of (6.8) is given as 



(6.9X 

Step Cvii): We use the transformation rule (2.11) to join the nonzero 
free angular-momentum lines j and J, 



C6.10) 
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By applying the transformation rule (2.12) to the zero angular-momentum 
lines in (6.10), we obtain the simple diagram 



((^ V')C ^ ^ 

^ + ^ To.t 7 

^ 5:^ 3- J 


( 6 . 11 ) 
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Step (vii): 





c. c ^ <r 


x[(2^+i:>{s.^^+0(x^-i-0'‘]~^(-f ^ ^ '^1 . 

H ^ 


C6.14) 
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Step (viii); By substituting (6,11) and C&*14l into (6.3), we obtain 
the antisymmetrized matrix element (the transition amplitude) for the 
C-AB;D Auger transition in a singly ionized closed-shell atom. 


Tv (c = r . J ' (15^ \ 2- 








3 ^c{ ^ ^ [s-ci 3-^ 3- ^ 


Cj5.15l 


Here we have used- the fact that and are all half-integers to 


simplify the phase factor. 
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APPENDIX A: ■ GLOSSARY OF TERMS 

Clebsch-Gordan coefficient, Wigner coefficient 

Alternative names for the vector-coupling coefficient. 

Wigner 3-j coefficient or symbol 

The symmetrized vector-coupling coefficient defined by Wigner (1951). 

Covariant 3-jm coefficient or symbol 

The vector-coupling coefficient defined in (2.2) in the covariant 
notation (Wigner, 1959). 

3n-j coefficient or symbol 

The 3-j symbol is defined as the triangular. delta in [2.19). The 
6-j and 9-j symbols have their usual meanings while the 12-j symbols 
and so on are not unique [Edmonds, 1957; Jucys et aT., 1962; El-Baz 
and Castel, 1972). 

3n-j coefficient or symbol of the first and second kind 

I 

The symmetrized recoupHng coefficient defined by Levinson and 
Vanagas (1957). Their definitions are given in [2.22). and [2.23). 
(Jucys et al., 1962; El-Baz and Castel, 1972). 

Subshel 1 

A collection of particle states having the same quantum numbers n, 
and j (or n and k). [See, e.g.. Grant, 1970.) 

Configuration 

A configuration is specified by the number of equivalent particles 
occupying each subshell and can be denoted by the aggregate {N }. 

‘ A 
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Antisymmetrized many-particle state 

The anti symmetrized state of a configuration. In this work, we mean 

, ^ *v 

exclusively the antisymmetrized jj-coupled state in which each subshell 
of equivalent particles is in a definite total subshell angular- 
momentum state and the whole many-particle system is in a definite 
grand total angular-momentum state. 

Anti symmetrized matrix element 

The matrix element between anti symmetrized many-particle states. 

Semi symmetrized many-particle state 

The many-particle state which. Is antisymmetric within each subshelT 
but is not antisymmetric with respect to exchange of two particles 
from different subshells. 

Semi symmetrized matrix element 

The matrix element between semi symmetrized -many-particle states. 

Active particles 

Those particles which actually participate in the interaction in a 
semi symmetrized matrix element. 

Spectator particles 

' Those particles which, as opposed to active particles, do not par- 
ticipate in the interaction in a semi symmetrized matrix element. 

jm-scheme matrix .element 

The. matrix element between uncoupled Dirac single-particle orbitals. 
Interaction graph 

^The geometric part of the jm-scheme matrix element of an interaction, . 
which depends only on the tensorial properties of the Interaction. 
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Interaction strength 

The dynamical part of the jm-scheme matrix element of an* interaction 
It can be expressed in terms of reduced matrix elements of tensor 
operators involved in the interaction. 

Bra and ket blocks 

The diagram blocks representing the angular-momentum coupling of the 
bra and ket states, respectively, of the many-particle system. 

Interaction block 

The diagram block representing the jm-scheme matrix element of an 
interaction, -including the interaction graph and the interaction 
strength. 

Spectator block 

The diagram block representing the scalar product of uncoupled 
spectator-particle states. 

Recoupling diagram 

The graphical representation of a recoupling coefficient. 
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APPENDIX B: GRAPHICAL FORMS OF OPERATORS AND POTENTIALS 

In this appendix we will consider the jm-scheme matrix elements of one- 
and two-particle operators, i.e., 

and 

(B.2) 

and their graphical representations. The Dirac orbitals in CB.l) and (B.2) 
are assumed to have the form 

II _ I / ' 

where the radial functions G and F are the large and small components, 

n< n< 

respectively, the are normalized two-component Dirac spinors. Here the 
orbitals are completely specified by the quantum numbers n, k, and m, which 
have their usual meanings (see, e.g.. Grant, 19.70). We will evaluate (B.l) 
and (B.2) in terms of radial integrals for various operators and potentials. 

We first define various functions, notations, and coefficients and 
present a few useful formulas. 


» 

(B.3) 
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(i) Different combinations of radial wavef unctions: 


W,b(r) - G,(r) Gb(r) + F,(r) Fb(r) , 

I'abf'') ” - fat'll ^b''') • 

''ab^’’> = '■’a<''> '"b^''> ^ fat'"' '^bt'"' ' 

“abt’’) = f’a''"' fbt'-) - fat'-) ■ 

. fabt'") = "abt'-) * 'abt'-l ■ 

W^)--“abt'")^4^'ab'^'- 

Cii) Various radial functions: 



where r^(r^) is the smaller (larger) of r^ and r^. 


CB.4) 


CB.5) 


CB,6) 


where and h^^ are the spherical Bessel and Hankel functions, respectively. 

I a) =e Cr, -rj ['R^^c r, r,:> r, ro] 

CB.7) 

with e(r-j - rg) the Heaviside step function ‘ 

r/ 

0 (z) = ■) 

^0 %<o , 
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rj 


^ ^ ^ = 




r > r 

i ^ ^JL > 


r. 




V 




^ijC^r,)-AjJwrJ, n<r^. 


(B.8) 


r 


jp-f 


oJ^r; 


-~F^ r, >Tj_ 

.' t -*'l ^ 


-|7- ip W r ) -^ C w r, ) , I", < rj. • 


I 


Cs.a), 
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(iil) The following notations are used to denote integrals: 


'^0 

rfG?^^(r,)R^Cr;rJ , 

^ ^ A 


/ V 

\}> Cr,) R.Cr, ri')Q<rJ) = 'Tf(4 j-iOTrfJJt 3- A) 

AC 3* M / 


,co /*^ 


K f ( ®lr;'p(rT)?(r,r,)Q (a),, 
/> A 3- M 


<? 




^ cn > = ir (4 ^ A ) ' -?0 


00 


Jo "^0 


where the parity selection function is defined as 


it(4a-4)= I 


4 Jl^ o4<< . 


0 
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(iv) Define the coefficients 

f 




(B.IO) 






7*4 


\ * % 


j- ^Z" i-i. i- 3; 


X 


>'4. ° -'^'^'L 0 -'4 


) 


(B.ll) 


c_M» S'"' ‘X!: “ X ^ s-iai 


^w\ 


Cv) 


3 I • C * K 4 ^ 


, Y/ n.^ 

KoJWva 4'’^ 


X)C • (^ ■ O • 


(B.13) 
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(vi) The vector-spherical -harmonics expansion of a U^:, 




(B.14) 


where 


(f>.- . = , 

Mi-0 ' 


(h = ''Ti \4.t^O 

s> 


Cu+0/2^^0} . 




CB.15) 


This was first derived by Mann and Johnson 0971), and we have verified it 
independently. Since no derivation has previously been given for. this ex- 
tremely useful formula, we present its derivation in Appendix C. Many of 
the techniques of the vector-spherical -harmonics expansion may be found in 
various books (Rose, 1957; Edmonds, 1957; Akhiezer and Berestetskii, 1965). 
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(vii) The vector-spherical -harmonics expansion -of B(l), 


6 CO =• V, tvl Ci.) O/c^^r 


is given as (Huang, 1978a) 




CB.16) 


where 




3+1 ‘'a 

t/'. . = 0 

■ 44 >4 




/aQ+O 






In the limit oj ^ 0, we have 




CB-17) 


and the radial function ip. in CB.17) becomes 

J * 
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oi) = n 




"pjo- 

“ I a 'I 




't ^(i+oV-^R- 

V^'+O ^ \--3't-3 ^1-1 W Z 3 W /^ 


(B.18) 


One-particle Operator 

We will present the jm-scheme matrix elements of the operator fCrl, 


-y ^ ^ 1 *J IC • r * i 

e» o*p, a, r, and a e in turn. 


(i) 


'<A|.fC0|k>=S^ . ■ 


ie^ 

— > 


CB.19) 


where the interaction strength, denoted by the cross "X“ in (B.19}, is 

X U; b) ’= <C ^lb ■ 


Note that ^ s 6, . s, 


'^b ’^a^b Vb ’ 


(8.2QJ 
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where 


<a|el b> = ■ 




CB.21) 


(ill) The matrix element of a*p may be obtained by making use of 
the formula 

U)= -i;<a.((x-r I b''> 


with the understanding that |b'> denotes the Dirac orbital with the 


substitutions 






. dr- ^ 


Tt - 

T,.Cr) = ( 

c\ 

- o(r 

J^) T, (r) 


for the radial parts of the orbital [b>. By using the expansion (B,14) and 
the formula 


we can easily obtain the result 


^(K C <Tu ,/N 


CB.26) 
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This agrees with the result obtained using other methods (see, e.g., Rose, 
1957; Grant, 1970). Hence the. graphical representation is 

f I ^ f 

^ (B.27) 


with the interaction strength 




CB.28) 


(iv) By using the expansion CB.14) and the formula CB-13), we obtain 


(B.29) 


where the radial integral is 


I 

1 


t e 


J2-»/€-+v 


<x.y 


Ca.3Q), 


Here e^ is the spherical unit vector CEdmonds, 1957). 
graphical representation is 


Therefore its 




4 


■ 4 “ 


Tl 


Ca.311 
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with the interaction strength 


t) ~ C^(0.; i>) (^l ; ip). 


(B.32) 


(v) For <a|r|fa>we obtain the same interaction diagram as (B.31) with 
the radial integral given by 




“ r > . 


CB.33) 


(vi) For <ala e^^*^lb> we use the expansion CB.14) and the familiar 
Rayleigh expansion of a. plane wave, 


t 


A 




1^0 






(B.34) 


The rest of the calculation is straightforward, and the result is 


ct- 








(B.35) 


Here the interaction graph is 




3a. 'f- 

f — > 


± 


n 




CB.36) 
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and the interaction strength is 

A 

% 




(B.37) 


‘ with D,. given as 
JJi 


3X 




. (3)''^(K..+ KbKV^, , 


. Otflcj 


'4 


' (j+O ^3'+i> , .^= a> 


CB.38) 


Two-particle Operator 


We win deal only with rotational invariant interactions, which are 
generally linear combinations of zero-rank tensors obtained by contracting 
tensors of the same rank. The matrix element of these interactions can be 
calculated with the result (Huang, 1971, 1978a, and 1978b} 




(B.39) 


where the interaction graph is 
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a ^ h\c, Sc ^ )rt. rt' -,■ J 


b 3^ 3 



^ .3, 


and the interaction strength is 


y.(«.b;c4) •= C-(o^i>\C(^:> X i<xh ^ cc( ), 
^ . J .3 


Here I.Cabjcd) is defined in terms of radial integrals, depending on the 
Specific form' of V(r-| 2 ). We summarize the results for various potentials: 

Ci) Coulomb potential: . 


X \ a ) ^ 

1 ^ ^^<KC a Lei 


(K.C .i t, <1 


^ 6 * 

(ii) Covariant photon interaction: 0 - — ~ — 

' ^ , *^12 
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- 0- C + V/,) 


odt^ 


£.\i&yi. 




(B.43) 


This Is obtained by using the expansion CB.14) and the formula 


1‘ 


'V-*> 


‘i 2^ ", ■' 'At; 

^+J^4w\+ I 




= (o 


f , s , r . 

^3- J?? J?A- 


(B.44) 


(ill) Transverse photon interaction: 

Cu»r 

. e ' 


- (c<,-0 "7“ I WV,^ J . 


(«.t:<La} = -Cl- VKcX Ct.+ '^<')fZ7T7'^ 


Oc>^ 


, • , . ^ . ShiZJr^ - 4iV42^ 

(B.45I 


wC,V«>\ 
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This is obtained by making use of the expansions CB.14) and CB.16) and the 
orthogonality relation 

(B.46) 






H-i Jit "ViVM 


(iv) Breit interaction: 


Tr, p ~J > 

^ ‘la. 


which represents the transverse photon interaction in the limit tu ^ 0. 


1 . (.b ; cM) f ^ ^77777) 

I £L\/AM I 

^ < S U/- 2)-p,^ > + -<p^^ u. ^ 




(B.47) 
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(v) Operator a-j'*a 2 » which is proportional to the leading imaginary 
part of the transverse photon interaction. 


a 


even- 






> 




(B.48) 
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APPENDIX C: VECTOR-SPHERICAL-HARHONICS EXPANSION OF u! 

From the properties of the two-component Dirac spinor 
prove that 


^ » r ~ XL > 


O' • L = - (,t<+0 j 


Hence three useful formulas may be derived, f.e.. 


( 


(A a. XL IZ 


3')v^ 3 


I' 






I. a 

d 


«^5 

we can easily 


Cc.l! 


• !<’«- >■ 


j<Am. 


t ^ } , 


LC.21 
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where etc., are defined in (B.13). To obtain the vector^ 

spherical -harmonics expansion of ut a Ul, i.e., 

a D 


= X Cr)Y cr)^ 


we simply evaluate the expansion coefficient as 


(C.3) 







(C.4) 


By using the formulars (C.2) and noting the relations 


Y -(j+i)’^c:(rxL )] Y. , 

A *»1 


4 


Y 


H3 ^ 


• '7"’4 - . 

= ^ T' ’ 

^ y A 


Qyyy 


K'Crxt)) 


and 


oi 


__ / 0 ^ \ 
~ V 9S O J ^ 


we can performs the angular Integration of CC-4). The expansion (B.14) of 
Og a is thus obtained. 
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